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Abstract 

We report results on studies of CP violation in the three-body charmless decay — > K^ir^ir^ . 
Evidence at the 3.9a level for large direct CP violation in B ± — > p(770)°K ± is found. This is the 
first evidence for CP violation in a charged meson decay. The analysis is performed using Dalitz 
analysis technique with a data sample that contains 386 million BB pairs collected near the T(4S) 
resonance, with the Belle detector operating at the KEKB asymmetric energy e + e~ collider. 

PACS numbers: 13.25.Hw, 11.30.Er, 14.40.Nd 
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INTRODUCTION 



Decays of B mesons to three-body charmless hadronic final states may provide new pos- 
sibilities for CP violation searches. In contrast to decays to two-body final states where 
direct CP violation can only manifest itself as difference in decay rates for B and B mesons 
to charge conjugate final states, in three-body decays it can also be observed as a differ- 
ence in relative phases between two quasi-two-body channels. A necessary condition for 
observing direct CP violation in a two-body decay is a non-trivial strong phase difference 
between the CP-conserving and CP-violating amplitudes contributing to a particular final 
state. Although this condition (if satisfied) also enhances the sensitivity to CP violation in 
three-body decays, it is not required in general and direct CP violation in quasi-two-body 
decays can also be observed with any strong phase difference via the interference with a 
nearby quasi-two-body or non-resonant amplitude(s). Although direct CP violation has 
been observed in decays of neutral K mesons |l| and recently in neutral B meson decays 
no CP violation in decays of charged mesons has been found to date. However, large direct 
CP violation is expected in some quasi-two-body modes Several other ideas to study 
CP violation utilizing decays to three-body final states have been proposed 0, 

First results on the amplitude analysis of the B^ — > K n ir^ decay are described in 
Refs. (1, the first results on searches for direct CP violation from independent fits of 



the B~ and B + samples are given in Ref. |10l. 111!]. The analysis of direct CP violation in 
the decay B ± — > K ± w ± Tt T described in this paper is based on a simultaneous fit to B~ and 
B + events. The results are obtained with a data sample of 357 fb -1 containing 386 million 
BB pairs, collected with the Belle detector operating at the KEKB asymmetric-energy e + e~ 
collider 13| with a center-of-mass (cm.) energy at the T(45) resonance (on-resonance data). 
The beam energies are 3.5 GeV for positrons and 8.0 GeV for electrons. For the study of the 
e + e _ — ► qq continuum background, we use data taken 60 MeV below the T(4S) resonance 
(off-resonance data). 



THE BELLE DETECTOR 

The Belle detector (l3j is a large-solid-angle magnetic spectrometer based on a 1.5 T 
superconducting solenoid magnet. Charged particle tracking is provided by a silicon vertex 
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detector and a 50-layer central drift chamber (CDC) that surround the interaction point. 
Two inner detector configurations were used. A 2.0 cm beampipe and a 3-layer silicon vertex 
detector was used for the first sample of 152 million BB pairs, while a 1.5 cm beampipe, a 4- 
layer silicon detector and a small-cell inner drift chamber were used to record the remaining 
234 million BB pairs ^]|- The charged particle acceptance covers laboratory polar angles 



between 9 = 17° and 150°, corresponding to about 92% of the total solid angle in the cm. 
frame. The momentum resolution is determined from cosmic rays and e + e~ — > fi + [i~ events 
to be (Jp t /pt = (0.30 © 0.19p t )%, where p t is the transverse momentum in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the CDC, an ar- 
ray of 1188 aerogel Cerenkov counters (ACC), and a barrel-like array of 128 time-of-flight 
scintillation counters (TOF); information from the three subdetectors is combined to form 
a single likelihood ratio, which is then used in kaon and pion selection. At large momenta 
(> 2.5 GeV/c) only the ACC and CDC are used to separate charged pions and kaons since 
here the TOF provides no additional discrimination. Electromagnetic showering particles 
are detected in an array of 8736 CsI(Tl) crystals (ECL) that covers the same solid angle as 
the charged particle tracking system. The energy resolution for electromagnetic showers is 
a E /E = (1.3 © 0.07/£ © 0.8/E 1 / 4 )%, where E is in GeV. Electron identification in Belle is 
based on a combination of dE/dx measurements in the CDC, the response of the ACC, and 
the position, shape and total energy deposition (i.e., E/p) of the shower detected in the ECL. 
The electron identification efficiency is greater than 92% for tracks with pi a b > 1.0 GeV/c 
and the hadron misidentification probability is below 0.3%. The magnetic field is returned 
via an iron yoke that is instrumented to detect muons and K\ mesons. We use a GEANT- 
based Monte Carlo (MC) simulation to model the response of the detector and determine 
its acceptance 



EVENT RECONSTRUCTION 

Charged tracks are selected with a set of track quality requirements based on the number 
of CDC hits and on the distances of closest approach to the interaction point. We also 
require that the track momenta transverse to the beam be greater than 0.1 GeV/c to reduce 
the low momentum combinatorial background. For charged kaon identification we impose a 
requirement on the particle identification variable which has 86% efficiency and a 7% fake 
rate from misidentified pions. Charged tracks that are positively identified as electrons or 
protons are excluded. Since the muon identification efficiency and fake rate vary significantly 
with the track momentum, we do not veto muons to avoid additional systematic errors. 

We identify B candidates using two variables: the energy difference AE = E B — E^ eSim = 

V ° 2 Pi + c4m ?) -^kam) and tne beam constrained mass M bc = ^\/E* 2 eam - c 2 (X^Pi) 2 > 
where the summation is over all particles from a B candidate; p« and rrii are their cm. three- 
momenta and masses, respectively. The signal AE shape is fitted by a sum of two Gaussian 
functions with a common mean. In fits to the experimental data, we fix the width and the 
relative fraction of the second Gaussian function from MC simulation. The common mean 
of the two Gaussian functions and the width of the main Gaussian are floated. The AE 
shape for the qq background is parametrized by a linear function. The Mb c distribution 
for the signal events is parametrized by a single Gaussian function. The M^ c width is 
about 3 MeV/c 2 and, in general, does not depend on the final state (unless photons are 
included in the reconstructed final state). The background shape is parametrized with an 
empirical function f(M^ c ) oc xy/l — x 2 exp[— £(1 — x 2 )] [l^, where x = Mio C /E^ e&m and £ is 
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a parameter. 



BACKGROUND SUPPRESSION 

The dominant background is due to e + e~ — > qq (q = u, d, s and c quarks) continuum 
events that have a cross-section about three times larger than that for the e + e~ — > T(4S) — > 
BB. This background is suppressed using variables that characterize the event topology. 
Since the two B mesons produced from an T(4S") decay are nearly at rest in the cm. frame, 
their decay products are uncorrelated and the event tends to be spherical. In contrast, 
hadrons from continuum qq events tend to exhibit a two-jet structure. We use ^hr, which 
is the angle between the thrust axis of the B candidate and that of the rest of the event, 
to discriminate between the two cases. The distribution of | cos 6W | is strongly peaked near 
Icos^thrl = 1-0 for qq events and is nearly flat for BB events. A Fisher discriminant is 
utilized for the further suppression of the continuum background. When combined, these 
two variables reject about 98% of the continuum background in the B ± — > K ± n ± it T decay 
while retaining 36% of the signal. A detailed description of the continuum suppression 
technique can be found in Ref. p and references therein. 

Another background originates from other B meson decays. We study the BB related 
background using a large sample of MC generated BB generic events. Note that charmless 
hadronic B decays that proceed via b — > s(d) penguin and b — > u tree transitions are not 
included in the generic MC sample and are generated separately. We find that the dominant 
BB related background to the K + it +r K~ final state is due to B + — > D°tt + , D° — > K + n~ and 
due to B + -> J/i)(i)(2S))K + , J/if>(i/j(2S)) -> decays. We_veto B+ -> D°ti + events by 

requiring \M(Kir) - M D \ > 100 MeV/c 2 . We also veto B + -> D°K + , D° -> tt+tt" signal by 
requiring \M(n + n^) — Md\ > 15 MeV/c 2 . To suppress the background due to tt/K misiden- 
tification, we also exclude candidates if the invariant mass of any pair of oppositely charged 
tracks from the B candidate is consistent with the D° — > K + n~ hypothesis within 25 MeV/c 2 
(~ 4a), regardless of the particle identification information. Modes with J /ip(ip(2S)) in 
the final state contribute due to muon-pion misidentification; the contribution from the 
J /ip(ip{2S)) — > e + e~ submode is found to be negligible after the electron veto requirement. 
We exclude J '/ip(ip(2S)) background by requiring \M(TT + n^) fl +^- — Mj/^\ > 70 MeV/c 2 and 
|M(7r + 7r~) M + M - — M^{2s)\ > 50 MeV/c 2 , with a muon mass assignment used here for the 
pion candidates. Yet another small but clearly visible background is due to B + — > J/ipK + , 
J/ip — > decay with a complicated series of particle misidentifications; the charged 

kaon from the B is misidentified as a pion, the n + is misidentified as a kaon and the yT 
as another pion. This background is excluded by applying a veto on the invariant mass of 
oppositely charged kaon and pion candidates: \M{K + it~) IJi + (l - — Mjm\ > 20 MeV/c 2 . The 
most significant background from charmless B decays is found to originate from B + — > rj'K + 
followed by rj' — > 7r + 7r _ 7. Another contribution comes from B + —> 7r + 7r + 7r _ decay, where one 
of the two same charge pions is misidentified as a kaon. Finally, we consider a background 
from the decay B° — > K + 7t~. Although it does not directly contribute to the AE signal 
region, this background should be taken into account in order to correctly estimate the qq 
component of the background. 
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FIG. 1: AE distributions for (a) B ± -> K^^if? combined; (b) B~ -> K - 7r - 7r + and (c) B + -> 
K + tt + tt~ candidate events. Points with error bars are data; the open histogram is the fit result; 
the hatched histograms are various background components. 



THREE-BODY SIGNAL YIELDS 

The AE distribution for —>■ K^tt^tt^ candidates that pass all the selection require- 
ments are shown in Fig. Q In the fit to the AE distribution we fix the shape and normaliza- 
tion of the charmless BB background components from the measured branching fractions 
and known number of produced BB events. For the BB generic component we fix only the 
shape and let the normalization float. The slope and normalization of the qq background 
component are free parameters. Results of the fit are shown in Fig. where different com- 
ponents of the background are shown separately for comparison. There is a large increase in 
the level of BB related background in the region AE < —0.15 GeV. This is mainly due to 
B — > Dn, D — > Kixix decay. This decay mode produces the same final state as the studied 
process plus one extra pion that is not included in the energy difference calculation. The 
semileptonic decays B — ► D^n, D — ► Kfiu^ also contribute due to muon-pion misidentifi- 
cation. The shape of these backgrounds is well described by MC simulation. Results of the 
AE fits are given in Table 

For the analysis of quasi-two-body intermediate states that contribute to the observed 
B^ 1 — > K ir iz^ three-body signal, we define the B signal and sideband regions as shown in 
Fig. |21 Defined in this way, the Mb c — AE sidebands are equivalent to the following side- 
bands in terms of the three-particle invariant mass M(K-ktc) and three-particle momentum 
P(Kirir): 

0.05 GeV/c 2 < \M(Kmr) - M B \ < 0.10 GeV/c 2 ; P(Ktttt) < 0.48 GeV/c. 



TABLE I: Results of the fits of the AE distributions. 



Final state 


G\ <J2 


Fraction of the 


N qq 


N B b 


Signal Yield 






MeV MeV 


main Gaussian 


(events) 


(events) 


(events) 


B- - 


-> K tt ir + 


15.6 ± 0.6 35.0 (fixed) 


0.84 (fixed) 


16932 ± 275 


8639 ± 226 


2248 ± 79 


B+ - 


-> K + ir + ir~ 


15.0 ± 0.6 35.0 (fixed) 


0.84 (fixed) 


17268 ± 274 


8828 ± 227 


2038 ± 76 


B±- 




15.3 ± 0.5 35.0 (fixed) 


0.84 (fixed) 


34188 ± 386 17452 ± 320 


4286 ± 99 
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FIG. 2: (a) Two-dimensional AE versus M^ c plot for the — ► iiT =l= 7r = ' = 7r =F signal MC events; 
(b) Definition of the signal and sideband regions on the M^ c — AE plane. 



and 



\M{Kim) - M B | < 0.10 GeV/c 2 ; 0.48 GeV/c < P^tttt) < 0.65 GeV/c. 



The signal region is defined as an ellipse around the M^ c and AE mean values: 



+ 



AE 2 



(M bc - M B ) 

(7.5 MeV/c 2 ) 2 ' (40 MeV) 



< 1, 



The efficiency of the requirements that define the signal region is 0.927; the total number 
of events in the signal region is 7757. The relative fraction of signal events in the B signal 
region is determined to be 0.512 ±0.012. There are 27855 events in the sideband region that 
is about seven times the estimated number of background events in the signal region. 



AMPLITUDE ANALYSIS 

The amplitude analysis of three-body B meson decay reported here is performed by means 
of an unbinned maximum likelihood fit. Details of the analysis technique are described in 
Ref. H. One of the important questions that arise in unbinned analysis is the estimation of 
the goodness-of-fit. As the unbinned maximum likelihood fitting method does not provide a 
direct way to estimate the quality of the fit, we need a measure to assess how well any given 
fit represents the data. To do so the following procedure is applied. We first subdivide the 
entire Dalitz plot into 1 (GeV/c 2 ) 2 xl (GeV/c 2 ) 2 bins. If the number of events in the bin is 
smaller than iV m i n = 25 it is combined with the adjacent bins until the number of events 
exceed the minimum required level. After completing this procedure, the entire Dalitz plot 
is divided into set of bins of varying size, and a x 2 variable for the multinomial distribution 
can be calculated as 
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M 2 (K + 7T) (GeV 2 /c 4 ) M 2 (K + 7i ) (GeV 2 /c 4 ) 



FIG. 3: Dalitz plot for the K^ir^ events in the (a) M^ c — A.E sideband region and in (b) B 
signal region. 

where is the number of events observed in i-th bin, and Pi is the number of predicted 
events from the fit. For a large number of events this formulation becomes equivalent to the 
usual one. Since we are minimizing the unbinned likelihood function, our "x 2 " variable does 
not asymptotically follow a \ 2 distribution but it is bounded by a \ 2 variable with (A^bms — 1) 
degrees of freedom and a x 2 variable with (iVbins —k — 1) degrees of freedom, where k is the 
number of fit parameters. Because it is bounded by two \ 2 variables, it should be a useful 
statistic for comparing the relative goodness of fits for different models. 



Fitting the Background Shape 

Before fitting the Dalitz plot for events in the signal region, we need to determine the 
distribution of background events. The background density function is determined from an 
unbinned likelihood fit to the events in the M^ C — AE sidebands defined in Fig. El Figure Efa) 
shows the Dalitz plot for sideband events. 

We use the following empirical parameterization to describe the distribution of back- 
ground events over the Dalitz plot 

Bk^s^Sk) = ttie"^ + a 2 e~^ + a 3 e~^ 

+ a 4 e- (/34Sl3+/3sS23) + a 5 e~ {/36Sl3+(37Sl2) + a 6 e- (/38S23+ft5Sl2) 

+ 7l \BW(K* (892))\ 2 + l2 \BW(p(770))\ 2 } (2) 

where s 13 = M 2 (K + ii~), s 2 3 = M 2 (n + n^) and «j («i = 1.0), $ and 7, are fit parameters; 
BW is a Breit-Wigner function. The first three terms in Eq. (j2J) are introduced to describe 
the background enhancement in the two-particle low invariant mass regions. This enhance- 
ment originates mainly from e + e~ —>■ qq continuum events. Due to the jet-like structure 
of this background, all three particles in a three-body combination have almost collinear 
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T 




FIG. 4: Results of the best fit to the i ; C ± 7r ± tt^ events in the M\, c — A.E sidebands shown as 
projections onto two-particle invariant mass variables. Points with error bars are data; histograms 
are fit results. 

momenta. Hence, the invariant mass of at least one pair of particles is in the low mass 
region. In addition, it is often the case that two high momentum particles are combined 
with a low momentum particle to form a B candidate. In this case there are two pairs with 
low invariant masses and one pair with high invariant mass. This results in even stronger 
enhancement of the background in the corners of the Dalitz plot. This is taken into ac- 
count by terms 4 — 6 in Eq. (J2j). To account for the contribution from real i\"*(892)° and 
p(770)° mesons, we introduce two more terms in Eq. (J2J), that are (non-interfering) squared 
Breit-Wigner amplitudes, with masses and widths fixed at world average values [17| . 

The projections of the data and fits for the background events are shown in Fig. The 
xV^bins value of the fit is 127.6/120. 



Fitting the B ± -> K ± vr ± 7r =F Signal 

The Dalitz plot for K ± k ± 'k t events in the signal region is shown in Fig. EJb). There 
are 7757 events in the signal region that satisfy all the selection requirements. As found in 
Ref. the B + — > K + -k + -k~ decay is well described by a matrix element that is a coherent 
sum of #*(892)°7r+ if *(1430)°7r+, p(770)°K + , f (980)K+, f x (13QQ)K + , Xc oK + quasi-two- 
body channels and a non-resonant amplitude. The fx (1300) K + channel is added in order to 
describe an excess of signal events at M(n + n^) ~ 1.3 GeV/c 2 . With current statistics, the 
contribution of f 2 (l270)K + is found to be significant, but not sufficient to fully explain the 
excess of signal events in this mass region. In this analysis we modify the model by adding 
two more quasi-two-body channels: f2{1270)K + and u(782)K + and change the parameteri- 
zation of the / (980) lineshape from a standard Breit- Wig ner function to a coupled channel 
Breit-Wigner also known as Flatte parametrization [18J. Although the u;(782) — > 7r + 7r~ 
branching fraction is only (1.70 ±0.27)% ^3], the u;(782) natural width is rather narrow. As 
a result a numerical factor of r(p(770)°)/r(u;(782)) ~ 18 is introduced in the B — ► uj(782)K 
amplitude (relative to the B — > p(770)K amplitude) that compensates the smallness of the 
d;(782) — > 7r + 7r~ branching. As the independently measured B ± — ► u(782)K ± branching 
fraction [19] is comparable to that for B ± — > p(770)° , the interference between these two 
amplitudes might significantly distort the p(770)° lineshape. Finally, for CP violation stud- 
ies the amplitude for each quasi-two-body channel is modified to include two components: 
one that is independent of the sign of the B charge and a component that changes sign with 
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FIG. 5: Results of the fit to events in the signal region. Points with error bars are data, 

the open histograms are the fit result and hatched histograms are the background components. 
Inset in (a) shows the K*(892) — Kq (1430) mass region in 20 MeV/c 2 bins. Inset in (b) shows the 
XcO mass region in 10 MeV/c 2 bins. 



the charge of the B meson. The resulting decay amplitude reads as 

M{B ± -> ir ± 7T ± 7T T ) = 

= a K *e iS «* (1 ± b K *e i<fi «*)A(K*{892)°) + a K *e iS ^ (1 ± b K */* K t )A{K*{U30)°) 
+ a p e iS »(l ± 6 p e^)^(p(770)°) + a^{l ± b u e iVu )A(u(7S2)) 
+ a fo e i5 fo(l ± 6 /o e^o)^ Flatte (/ (980)) + a h e i5 ^(l ± 6 A e^)^(/ 2 (1270)) 
+ a fx e a 'x (1 ± b fx e*"*)A{fx) + a^e"** (1 ± b Xc0 e^)A( Xc o) 

+ Ar(^ ± VT ± 7T T ) (3) 

with the non-resonant amplitude A nT parametrized by an empirical function 

Ar(^ ± vr ± 7r T ) = aYe- asi3 e iS ™ + a^e"^ 23 ^ , (4) 

where S13 = M 2 (K ± 7r T ), S23 = M 2 (_7rJ~7r _ ). Note that alternative parameterizations of the 
non- resonant amplitude possible [1, [ll|. The amplitudes and 6«, relative phases 5i and ipi, 
mass, g wn and qkk of the / (980), mass and width of the /x(1300), and parameter a of the 
non-resonant amplitude are fit parameters. With such a parameterization of the amplitude, 
the CP violating asymmetry Acp for a particular quasi-two-body channel can be calculated 
as 

_ B{B~ ->/-)- -» /+) _ 2b cos y> 

CPU" J ^ J-) ^ J+) I -/ 

To reduce the number of free fit parameters, we fit the data in two steps. First we 
fix all 6j = and fit the data assuming no CP violation. From this fit we determine 
the parameters of the /y(1300) (M(/ x (1300)) = 1.449 ± 0.013 GeV/c 2 , r(/ x (1300)) = 
0.126 ± 0.025 GeV/c 2 ), / (980) (M(/ o (980) = 0.950 ± 0.009 GeV/c 2 , = 0.23 ± 0.05, 
gKK — 0.73 ± 0.30) and the parameter of the non-resonant amplitude (a = 0.195 ± 0.018). 
We then fix these parameters and repeat the fit to data with bi and ipi floating. In addition, 
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FIG. 6: Helicity angle distributions for i^ ± 7r ± 7r =F events in different regions: (a) K*(892)° 
(0.82 GeV/c 2 < M{K + tt~) < 0.97 GeV/c 2 ); (b) #3(1430) (1.00 GeV/c 2 < 
M(K+tt~) < 1.76 GeV/c 2 ); (c) p°(770) (0.60 GeV/c 2 < M(^+^-) < 0.90 GeV/c 2 ) and 
(d) / (980) (0.90 GeV/c 2 < M(7r + 7r _ ) < 1.06 GeV/c 2 ). Points with error bars are data, the open 
histogram is the fit result and the hatched histogram is the background component. 



we also assume no CP violation in B ± -> u^l%2)K ± and for the non-resonant amplitude. 
Possible effects of these assumptions are studied and considered in the final results as a part 
of the model uncertainty. 

The numerical values of the fit parameters are given in Table ITT1 The xV^bins value of 
the fit is 182.5/141 with k — 32 fit parameters. Fit projections and the data are shown 
in Fig. El Figure El shows helicity angle distributions for several regions, where the helicity 
angle is defined as the angle between the direction of flight of the ir~ in the h + ir~ rest frame 
and the direction of B + candidate in the h + n~ rest frame. Gaps visible in Fig. El are due to 
vetoes applied on invariant masses of two-particle combinations. All plots shown in Figs. El 
and El demonstrate good agreement between data and the fit. 

The statistical significance of the asymmetry quoted in Table |H] is calculated as 
a/— 2 ln(£o / £max) , where £ max and Cq denote the maximum likelihood with nominal fit 
and with the asymmetry fixed at zero, respectively. The only channel where the statistical 
significance of the CP asymmetry exceeds the 3cr level is — > p(770)°ii' ± . Figures |7(a,b) 
show the 7r + 7r~ invariant mass distributions for the p(770)° — /o(980) mass region separately 
for B~ and B + events. The effect is even more apparent when the M(tt + 7t~) distribution 
for the two helicity angle regions (cos9]f < and cos9]f > 0) shown in Fig. Efc-f) are 
compared. 
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TABLE II: Results of the best fit to if ± 7r ± 7T F events in the B signal region. The first quoted error is statistical and the second is the model 
dependent uncertainty. The quoted significance is statistical only. 



Channel 


CP averaged 


5, 


b 




A C p, 


Significance, 




fraction, % 


degrees 




degrees 


% 


a 


K*(892)°7T ± 


13.0 ± 0.8lg;f 


(fixed) 


0.078 ± 0.0331^03 


-18 ±441^3 


-14.9±6.4l°| 


2.6 


ifo(1430) 7r ± 


65.5 ± 1-5±§;| 


55 ± 4+1 


0.069 ± 0.03llg;gjg 


-123±16lJ 


+7.5 ± 3.8l^ 


2.7 




7.85 ± 0.931^9 


-21 ± U+\t 


0.28 ±0.111°;°* 


—125 ± 32±Jg 


+30 ± III4 1 


3.9 


^(782)^+ 


0.15 ± 0.12±g;g| 


100 ± 


(fixed) 








/o(980)K± 


17.7 ± 1.6+^ 


67 ±lllJ5 


0.30±0.19l° 1 ;°j 


-82 ± 8l^ 


-7.7 ±6.5+^ 


1.6 


/2(1270)A-± 


1.52 ± 0.35±g;§| 


140 ± ll +1 7 8 


0.37 ± 0.17±g;J| 


-24 ±291^ 


-59 ± 22+j 


2.7 


/x(1300)k+- 


4.14±0.81±g;|S 


-141 ± 10±| 


0.12 ±0.17l°$ 


-77 ± 56±H 


-5.4 ± 16.5l^° 4 3 


1.0 


Non-Res. 


34.0 ± 2.2+^ 


Sf = —11 ± 5±| 
5f = 185 ± 20l^ 


(fixed) 








Xc0^ ± 


1.12 ±0.12l°;g£ 


-118±24l^ 


0.15 ±0.35+^ 


-77±94ljf 


-6.5 ± 19.6+jj^ 


0.7 




M(Jt + 7i ) (GeV/c 2 ) M(Jt + 7i ) (GeV/c 2 ) 
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FIG. 7: M(tt +/ k~) mass spectra for B~ (left column) and B + (right column) for different helicity 
regions: (a,b) no helicity cuts; (c,d) cos 9jj < 0; (e,f) cos 8™ > 0; Points with error bars are data, 
the open histogram is the fit result and the hatched histogram is the background component. 



SYSTEMATIC & MODEL UNCERTAINTIES 

The dominant sources of systematic error are listed in Table ITTTl The systematic uncer- 
tainty in charged track reconstruction is estimated using partially reconstructed D* — > Dtt 
events and from comparison of the ratio of i] — > 7r + 7r~7r° to r] — » 77 events in data and MC. 
The uncertainty from the particle identification efficiency is estimated using pure samples of 
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TABLE III: List of systematic errors (in percent) for the three-body — ► K^tt^tt^ branching 
fraction. 



Source 


Error 


Charged track reconstruction 


3.0 


PID 


4.5 


Event Shape requirements 


2.5 


Signal yield extraction 


3.9 


Model 


1.2 


MC statistics 


1.0 


Luminosity measurement 


1.0 


Total 


7.4 



kaons and pions from D — > K~tt + decays, where the D flavor is tagged using D* + — > D ty + 
decays. The systematic uncertainty due to requirements on event shape variables is esti- 
mated from a comparison of the | cos# t hr| and T distributions for signal MC events and 
B + — > D°7T + events in the data. We estimate the uncertainty due to the signal AE shape 
parameterization by varying the parameters of the fitting function within their errors. The 
uncertainty in the background parameterization is estimated by varying the relative frac- 
tion of the BB background component and the slope of the qq background function within 
their errors. Reconstruction efficiency is determined using MC events distributed over phase 
space according to the matrix element corresponding to the best fit to data. The relevant 
systematic uncertainty is estimated to be at the level of one percent. Finally, to account for 
variations in reconstruction efficiency due to modifications in the detector setup and due to 
non- uniform data-taking conditions (mainly beam related background conditions), we gen- 
erate signal events with background events embedded. Background events are recorded with 
random triggers for each experiment. Signal MC events are generated for each experiment 
with statistics proportional to experimental data. The overall systematic uncertainty for 
the three-body branching fraction is 7.4%. 

Note that in asymmetry calculation most of these systematic uncertainties cancel. The 
few remaining sources are listed in Table ITVl Systematic uncertainty due to possible asym- 
metry in background from charmless B decays is estimated by introducing an asymmetry 



TABLE IV: List of systematic errors for CP violating asymmetry. 

Source SAcp 
Rare background 

p ^ +0.003/ - 0.004 

rfK± +0.001/ - 0.001 

K^ic* +0.004/ - 0.004 

Detector asymmetry +0.023/ — 0.023 

Signal yield extraction +0.011/ - 0.011 

Total +0.031/ - 0.029 
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equal to experimentally measured central value lj| increased by one standard deviation 
to each charmless background component one by one and refitting the data. The possible 
bias due to intrinsic detector asymmetry in reconstruction of tracks of different charges is 
estimated using B — > Dtt events in data. 

To estimate the model dependent uncertainty in the branching fractions and asymmetries 
for individual quasi-two-body intermediate states, we vary the default model and repeat 
the fit to data. Namely, we add one additional quasi-two-body channel which is either 
fr(1410)°7r + , ir(1680)°7r + , or fC;(1430) 7r + or remove ^(782)^+ or f 2 (1270)K + channel 
from the default model, use different assumptions on the spin of the /x(1300) state and 
use different parameterizations of the non-resonant amplitude. For estimation of the model 
uncertainty in charge asymmetries for individual quasi-two-body channels in addition to 
model variations we fit the data with different assumptions on CP violation in different 
channels. Finally, we check the consistency of the Acp results with those obtained from 
independent fits of B~ and B + subsamples. 



BRANCHING FRACTION & CHARGE ASYMMETRY RESULTS 



In the preceding section we determined the relative fractions of various quasi-two-body 
intermediate states in the three-body i? ± — ► K^n^rt^ decay. To translate those numbers 
into absolute branching fractions, we first need to determine the branching fraction for the 
three body decay. To determine the reconstruction efficiency, we use MC simulation where 
events are distributed over the phase space according to the matrix elements obtained from 
the best fit to data. The corresponding reconstruction efficiency is (22.4 ±0.2)%. Results of 
the branching fraction and CP-violating asymmetry calculations are summarized in Table IVl 



TABLE V: Summary of branching fraction results. The first quoted error is statistical, the second 
is systematic and the third is the model uncertainty. 



Mode 


B(B + - 


-+ Rh+ -> K+tt+tt-) 


B(B+ -> Rh+) 


A C p, 






xlO 6 


xlO 6 


% 



K ± 7r ± 7rT Charmless - 48.8 ±1.1 ±3.6 4.9 ± 2.6 ± 3.0 

^*(892)°[K+^-]7r+ 6.45 ± 0.43 ± 0.48±g;|| 9.67 ± 0.64 ± 0.72±g || -14.9 ± 6.4 ± 3.0±g;| 

K^(U30)[K + tt-}tt+ 32.0 ± 1.0 ± 2.4±J;J 51.6 ± 1.7 ± 6.8j£? ±7.6 ± 3.8 ± 3.0±g:g 

p(770)°[ir + TT-]K + 3.89 ± 0.47 ± 0.29^^9 3.89 ± 0.47 ± 0.29±g jg ±30 ± 11 ± S.O^ 1 

/o(980)[7r+7r-]iir + 8.78 ±0.82 ±0.651?;^ - -7.7 ± 6.5 ± 3.0^ 

/2(1270)[vr + 7r-]K+ 0.75 ± 0.17 ± 0.06±g;^ 1.78 ± 0.41 ± 0.14±g;^ -59 ±22 ±3.0^ 

Non-resonant - 16.9 ± 1.3 ± l.S+ll 

Xco[tt + tt~]K+ 0.56 ± 0.06 ± 0.04i°:o4 112±12±18^ 4 -6.5 ± 19.6 ± 3.0l^ 



DISCUSSION & CONCLUSION 

With a 357 fb _1 data sample collected with the Belle detector, we made the first analysis 
of direct CP violation in the three-body charmless decay B^ — > i ; T ± 7r ± 7r =F . Results on 
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branching fraction and CP- violating asymmetry calculations are summarized in Table 
In all except the — > p(770)°A' ± channel the measured asymmetry is below 3a statistical 
significance. Evidence for large direct CP violation the decay — > p(770)°iC i: is found in 
agreement with our results obtained with 253 fb _1 [HI and with results by BaBar 11 1. This 
is also in agreement with some theoretical predictions . The statistical significance of the 
asymmetry observed in B ± — > p(770)°K ± is 3.9a. Depending on the model used to fit the 
data the significance varies from 3.7a to 4.0a. If confirmed with a larger data sample, this 
would be the first observation of CP violation in the decay of a charged meson. 
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